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EMG Pattern Analysis of Isometric Muscular Contraction in Thermoregulatory Heat Pro- 
duction* 

(I) Cluster Analysis and Factor Analysis of Grouping Discharge 
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Hirokazu Fujimoto**, and Mitsuo Kosaka, (Department of Epidemiology and Environ- 
mental Physiology, Institute for Tropical Medicine, Nagasaki University) 


Abstract: In a neutral temperature environment, six healthy human volunteers aged 22- 
33 years were examined by computer analysis of EMG patterns of dominant biceps bra- 
chii muscle throughout its isometric contraction against a constant load. The isometric 
contraction was performed from a seated position with the elbow resting on a flat surface 
and 45° angle to the table (90° flexion). 

Fatigue EMG of unilateral biceps brachii was recorded during 40% maximum contrac- 
tion value (MCV) isometric contraction, using both surface and fine needle electrodes. 
Properties of grouping discharge in EMG patterns of cold shivering and muscular fatigue 
were investigated. The interval of consecutive two pulses, mean frequency and mean 
duration of grouping discharge were analyzed by Cluster analysis as well as by Factor 
analysis using an analogue computer ATAC 450 and a personal computer PC-8001. 

Two different modes of interval of grouping discharge were found from a dendrogram 
of P-mode Cluster analysis. One had an interval time of 1 一 19 msec and the other had 
20—50 msec of consecutive two pulses of grouping discharge. The data of Factor ana- 
ysis using a Direct Varimax method also revealed a significant difference between the 
two modes interval, namely in the one, less than 19 msec and in the other, more than 
21 msec, respectively. 

These results indicate that the range of interval of consecutive two pulses of grouping 
discharge in EMG of cold shivering and isometric muscular contraction is conclusively 
20 msec and that the range is identical with the result of previous studies. The dif- 
ference between Cluster analysis and Factor analysis was further discussed to evaluate 
the purposes and nature of the statistical analysis. 
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< は C め に > 


筋 収縮 こよ る 熱 産 生 機 構 の 一 つ で ある 寒冷 ふる え 
の 筋 電 図 波形 に みる 群 化 放電 19~20 (Grouping Dis- 
charge 以下 G.D. と 略す ) の 平均 周波 数 は 小 動 
物 ほ ど 大 きぐ, 同一 動物 種 で は , G. D. の 平均 周波 
数 と 体重 の 間 に は , 相関 の な いこ と が 報告 され て い 
2919, 先 に 我々 は 正常 タウ サギ の 寒冷 あぁ る を 筋 電 図 
上 の G.D. の 平均 周波 数 お よび 各個 の G.D. の 
平均 持続 時 間 を デー タ 処 理 シ ステ ム に て 解析 し た 結 
R, 連続 する 2 個 の 筋 放電 スペ イク 間隔 が 20msec 
以下 を G.D. の 条件 と する と , G.D. の 平均 周波 数 
は 22.5c/s と な り 従 来 の 報告 結果 19 と 有意 差 の 
な いこ と が わか っ た ゆめ . LEAT, ある え 筋 電 図 上 


に 出現 する G.D. は 筋 の 疲労 状態 の 筋 電 図 に も 出 


現す る こと が 知ら れ て いる 1 の そこ で 今回 等 尺 
性 筋 収縮 こよ る ヒト の 筋 の 疲労 状態 の 筋 電 図 に 観察 
され た G. D. を 同様 の デー タ 処 理 シ レス テム に て 解 
析 す る こと を 試み た . FOB, G.D. の 解析 の 前 段 
階 と し て , G. D. の 条件 設定 すなわち 2 つの 筋 放電 
スペ イク の 間隔 の 設定 を 多 変 量 解析 の 一 つ で ある , 
クラ スタ ー 分 析 ?52! の お よび 因子 分 析 ?7? を 用 いて 行 
O, その 結果 を 従来 の 報告 結果 "と 比較 検討 し , 
この 分 析 方 法 ? に つい て 幾ら か の 考察 を 加え た . 


<A 法 > 

(DERA 

測定 の 対象 と し た 被験 者 は , 6 人 の 健康 な 成人 男 
FC, その 年 令 構 成 は 22~33 才 (平均 28.3 オ ) で あ 
っ た . 
(② 実 験 方 法 

室温 25°C, 湿度 65 の 環境 条件 下 で 被験 者 を NL 
に 正 対 させ , 机 の 面 に 対し , 上 腕骨 軸 が 45?, 更に 
肘 を 90? 屈曲 させ た 状態 を と ら せ , 上 腕 二 頭 筋 の 等 
尺 性 収縮 を 行わ せしめ た Fig. 1 参照 ). 実験 に 先 
立ち , 各 人 の 上 腕 二 頭 筋 の 最大 随意 収縮 力 を バネ 計 
り に て 測定 し , その 40 を 負荷 量 と し た . 

実験 は , これ 以上 肘 を 90° に 保持 で き な く な っ た 
時 点 で 終了 と し た . 
(3 反応 指標 記録 

筋 電 図 は , 上 腕 二 頭 筋 よ り , 双極 針 電極 を 用 い 測 
EL, 増幅 し た 筋 電 図 波形 を シン クロ ・ ス ュー プ で 
観察 する と 共に , データ 処理 シス テム ATAC 450 


Fig. 1. Method for isometric contraction. A 
seated position with the elbow resting 
on a flat surface and 45° angle to the 
table (90° flexion). 


(日 本 光電 工業 KK), お よび パー ソナ ルコ ン ピ ュ 
ー タ PC-8001 (日 本 電気 KK) で 分 析 す る 目的 で 
デー タレ ュー ダ へ 収録 し た . 


E R> 


MRA (64) か ら 記 録 し た 筋 電 図 に つい て , & 
変量 解析 を 行っ た が , 以下 , その 典型 的 な 2 例 に つ 
いて 報告 する . 

A. 被験 者 [ (Cluster 分 析 に よる 分 類 ) 

① 収 録 さ れ た 筋 電 図 に つい て , 連続 する 2 つの 筋 
放電 パル ス の 間隔 が 1 て 50msec の イン ター バル を 
l msec 毎 に デー タ 処 理 シ ステ ム (ATAC 450) を 
用 い カ ウン ト し た . カウ ント は , 実験 開始 時 か ら 終 
了 時 まで , 1 分 毎 に 6 区 間 行 われ た . Fig. 2 は 最大 


NUN TM S 


e Ae pes 


20msec 


Fig. 2. EMG pattern changes during muscular 

fatigue 
EMG of localized muscular fatigue of 
unilaterel biceps brachii during a 4025 
maximum contraction value isometric 
contraction. 

A: Initial stage of fatigue 

B: Middle stage of fatigue 

C: Terminal stage of fatigue 


count nnmbers (100 count//DIV) 


1 min after 4 min after 


2 min after 5 min after 


3 min after 6 min after 


count numbers (100 count//DIV) 
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Fig. 3. Histograms of interval changes in consecutive two pulses of grouping 
discharg. EMG recorded in the time course of muscular fatigue of 
dominant biceps brachii during isometric contraction. 
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筋 収縮 力 の 40 多 負荷 に お ける 等 尺 性 筋 収縮 の A…… 
初期 1 分 経過 ), B… BER (3 分 経過 ),C…… 
終期 (6 分 経過 ) の 筋 電 図 , 及び Fig. 3 は , 2 つ 
の 筋 放電 パル ス 間 の イン ター バル の カウ ント 数 に よ 
り 作成 され た 1 分 毎 6 区 間 の ヒス ト グ ラ ム で ある . 
②2 つの 筋 放電 パルス 間 の 各 イ ンター バル の , E 
験 開始 か ら 終了 まで の カウ ント 数 の 1 分 毎 の 変化 を 
も と に , 1~50msec の 各 イ ンタ ー バ ル に つい て の 
相関 行列 を 求め , P-mode Cluster 分 析 を 行っ た . 
P-mode Cluster 分 析 の 結果 , Fig. 4 に 示す dendro- 
gram が 得 ら れ た . 

Fig.4 の dendrogram の 横 方 向 は 相関 係数 を 表 
わし , 右 へ 行く ほど 正 の 相関 が 高く な っ て いる , ま 


た 縦 方 向 は 計算 に より 分 類 さ れ た 2 つの 筋 放電 ペル 
ス 間 の 各 イ ンタ ー ベ ル を 表わし て いる . Similarity 
level of clustering (相関 係数 ) 一 0.412 で 分 類 し た 
所 , 各 イ ンタ ー ベ バル の 順序 に 若干 の 変動 が 見 られ る 
も の の , 1~19msec £20~50msec の 2 つの Clus- 
ter に 分 ける こと が 出来 た . この 2 つの Cluster の 
各 イ ンタ ー バ ル に つい て は , Cluster の 内 部 で は , 

各々 が 相対 的 に 正 の 相関 を な し , Cluster 間 に お い 
て は , 各々 が 負 の 相関 を な し て いる と 言え る . つま 
り , 2 つの パル ス 間 の イン ター ベル が 1 19msec 
の も の は , 筋 が 疲労 に 致 る まで の 頻度 の 時 間 に よ る 
ダイ ナミ ッ ク な 変化 が , 20—50msec の も の と は 逆 
の 動き を し て いる こと が 示唆 され る . 故に , G.D. 


Correlation coefficient 


—0.50 0.00 


+0.50 1.00 

1 msec 0.996 
2 0.998 
3 0.990 
4 0.999 
5 0.906 
6 0.964 
7 0.730 
8 0.996 
9 0.999 
10 0.986 
14 0.967 
11 0.980 
13 0.949 
12 —0.060 
15 0.704 
16 0.861 
17 0.482 
18 0.693 
19 —0.412 
20 0.987 
25 0.898 
22 0.960 
42 0.979 
49 0.816 
21 0.942 
23 0.983 
24 0.876 
26 0.993 
37 0.973 
30 0.955 
28 0.983 
33 0.747 
27 0.967 
50 0.936 
29 0.967 
41 0.925 
32 0.980 
44 0.900 
43 0.785 
31 0.978 
38 0.990 
40 0.945 
35 0.988 
39 0.967 
36 0.981 
48 0.905 
34 0.909 
46 0.880 
45 0.483 
47msec 


Fig. 4. Dendrogram of a subject | . Clusters divided into two parts at a 
correlation coefficient —0.412. 


の 条件 設定 値 . つ まり 2 つ の 筋 放 電 ぶ ルス 問 の イン 
ター バル が その 値 以 内 だ と 1 つの グル ー プ 内 と 見 な 
され , それ 以上 だ と 次 の グル ー プ に 移っ た と 見 な さ 
れる 値 を この 被験 者 に つい て は 20msec と 設定 し 
た だ た: 

B. RRA] (因子 分 析 に よる 分 類 ) 

DORRE | の 場合 と 同様 に , 実験 開始 か ら 終 了 ま 
で の 筋 電 図 に つい て , 2 つの パル ス 間 の 各 イ ンタ ー 
バル を 1 分 毎 に カウ ント し た . 

② 被 験 者 T の 場合 と 同様 に , 各 イ ンタ ー バ ル に つ 
いて の 相関 行列 を 求め , P-mode Cluster 分 析 を 行 
っ た . 計算 に ょ り 得 られ た Fig.5 の dendrogram 
に 示さ れ て いる よう に , 各 イ ンタ ー ベ ル が 3 以上 の 
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Cluster を な し , dendrogram を 見 た だ だけ で は 2 つ 
の グル ー プ 間 の イン ター ズバ ル を 設定 する 手がかり が 
得 ら れ な か っ た . この 理由 と し て は , この 被験 者 の 
2 つの 筋 放電 パル ス 間 の 各 イ ンター バル の 時 系 列 変 
化 に 複数 の 変動 要因 が 働い て いて , 各 イ ンタ ー バ ベル 
の 相対 的 関係 を 視覚 的 に 示す Cluster 分 析 で は , そ 
の 計算 結果 が 複雑 化し た こと が 考え られ る ・ 
③Cluster 分 本 を 行っ た だ け で は , 2 つの グルー 
プ 間 の イン ター バル を 設定 する 手がかり が 得 ら れ な 
か っ た の で , Cluster 分 析 で 用 いた , 2 つの 筋 放 電 
バル ス 間 の 各 イ ンタ ー バ ル に つい て の , 同じ 相関 行 
列 を 用 い , Direct Varimax 法 に よる 因子 分 析 を 行 
っ た . Table 1 は 計算 に より 得 ら れ た 因子 行列 で あ 


Correlation coefficient 


—0.50 0.00 


+0.50 1.00 

1 msec 0.861 
35 0.919 
40 0.648 
28 0.901 
32 0.959 
37 0.800 
36 0.517 
19 0.967 
30 0.920 
38 0.812 
24 0.921 
34 0.815 
33 0.192 
27 0.397 
29 0.777 
Al 0.996 
46 0.588 
43 0.969 
48 —0.137 
2 0.939 
14 0.536 
3 0.971 
5 0.929 
4 0.986 
6 0.758 
7 0.955 
8 0.892 
15 0.766 
9 0.974 
12 0.404 
45 —0.039 
10 0.900 
49 0.755 
11 0.672 
16 0.852 
42 — 0.202 
13 0.959 
47 0.501 
18 0.141 
17 0.877 
39 0.683 
20 0.859 
44 0.540 
21 0.854 
23 0.985 
25 0.811 
26 0.954 
31 0.367 
22 0.909 
50msec 


Fig. 5. Dendrogram of a subject | . Clusters not to be divided into two 
parts at any correlation coefficient 
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Table 1. Factor loadings by Direct Varimax 


method (for details see text). 


Intervals 
(msec) 


UNNA 64 Me 


10 


^ 


1 tA YA LA し 』 ted ty l 
ONO cn nk 


» 


LE 


Contributions 


Factor-loadings 


] 
—-0.5777 
+0. 6313 
+0. 8084 
+0, 9221 
+0,9076 
+0. 97554 
+0., 9933 
+0. 7808 
+0,8195 
+0, 46392 
+0, 2277 
+0, 8293 
~O.1851 
+0. 6742 
*0,.8928 
-6.0276 
+0. 1464 
-O, 0381 
-0, 97412 
*0.0140 
-0,4523 
-O,3179 
-9, 4806 
-9,9212 
-O, 3684 
ーQ 6591 
-0.18975 
-0.6910 
-O.6000 
-O,8733 
-0.63545 
—-0,49907 
-0.9277 
-O, 8502 
-0.5714 
-0.5078 
-0, 6288 
ーQ 。 9408 
+0,2726 
-O, 3208 
-O, 7865 
+0,0155 
+0. 125i 
+0, 2390 
+0. 5848 
-0.7750 
+0,0572 
-0,.10461 
+0, 3931 
-0.3119 


+0. 3950 


I 
+Ü. 4024 
+O. 4567 
-0.1225 
-0.2172 
*O0.0311 


^" m^ 
-0 a dle 


-0.0214 
*0.1626 
-O, 2586 
ーQ 。 6891 
-O, 5516 
-9, 0561 
+0,1398 
40.1870 
*0.2811 
—-0,.9913 
40,9351 
*0.9177 
-9, 0343 
+0, 402 

+0. 3777 
+4+0,1 313 
+0., 1448 
+O., 2106 
+0, 2545 
+0, 3228 
-0, 6354 
+0. 1859 
-O, 14235 
-0.1152 
+0, 6839 
-O,0881 
ー0 。 1537 
*0.497X 
+0. 0885 
+0., 5194 
+0. 1239 
0. 0644 
+0, 7399 
+0, 2860 
-O, 5022 
—0.9116 
—O0. 5103 
+0, 2574 
-9, 0589 
-O, 4554 
+9, 2891 
-9,4216 
-0.7202 


-0.2427 


*0.1910 


5. 
因子 行列 と は , 想定 され る いく つか の 因子 と 項 
BE, つま り 2 つ の 筋 放電 パル ス 問 の 各 イ ンタ ー ベ バル 
と の 相関 行列 で あり , 1~50 行 目 が 各 イン ター バ 
^, 1, 2 列 目 が , 第 1 因子 , 第 2 因子 で ある . こ 
こ で 因子 と 言う の は , デー タ の 変化 を 規定 する 変動 
要因 で ある が , その 実体 は , 各 項 目 と の 相関 係数 か 
ら 推 測 さ れる . 本 研究 の 場合 , 各 イ ンタ ー ズ ベル の 分 
類 が 目的 で ある の で , 因子 の 解釈 は , 特に 行わ な か 
っ た . 計算 で は , 第 2 因子 まで 求め て いる が , HA 
率 に つい て 見 て みる と , 第 1 因子 が 39.5 に 対し , 
第 2 因子 も 19.1 有 と か な り 高 い 値 を 示し て いる . 帰 
寄 率 と は , 想定 され た 因子 が , 変動 要因 全体 の 中 で 
どれ だ け 説 明 す る こと が 出来 る か と 言う 確率 の こと 
で ある が , 第 2 因子 の 帰 寄 率 の 高 さ が , Cluster 分 
析 の 計算 結果 を 複雑 化 さ せ た 原 因 で ある こと が 示唆 
され る . 

第 1 因子 に つい て 見 て みる と 1 msec~10msec ま 
で と 12, 14, 15msec に 正 の 相関 が 見 られ , 11, 13, 
16~18, 20msec は 相関 が 見 られ ず , 19msec と 21 
msec 以上 に は , 一 部 例外 が ある も の の 負 の 相関 が 
見 られ た . Lot, 負 の 相関 と 無相 関 の 変 化す る 20 
msec を , この 被験 者 の 。 2 つの グル ー プ 間 の イン 
ター バル に 設定 し た . 


以上 , 被験 者 1・ を 用 いた 誘発 等 尺 性 筋 収 縮 筋 
電 図 波形 上 に 観察 され る 筋 放 電 パ ルス の 間隔 は 20 
msec 区 特異 性 が 認め られ る 事 が Cluster 分 析 お よ 
び 因 子 分 析 に よっ て 証明 され た . 


CX > 


筋 電 図 上 の G.D. の 分 析 に お いて は , grouping 
の 条件 つま り 2 つの グル ー プ 間 の イン ター バル の 設 
定 が 重要 な 問題 と な た っ て くる が , 従来 の 研究 *-1 
19770 に お いて は , その 判断 に 主観 的 要素 が 入り 込 
む 可 能 性 が 大 きか っ た . これ に 反し 条件 設定 に 多 変 
量 解 析 の 結果 を 手がかり に する と , 判断 が , より 客 
観 的 と な り , 本 研究 に お いて 有効 な 手段 で ある と 考 
えら れる . 中 で も Cluster 分 析 は , 各 項 目 つ まり , 
2 つの パル ス 間 の イン ター バル 相 互 の 相対 的 関係 を 
dendrogram と 言う 形 で 視覚 的 に 示す こと が 可能 で 
あり , また 因子 分 析 は , 各 項 目 と , 想定 され る 因子 
と の 関係 を , 相関 係数 で 示し , 変動 要因 が 2 つ 以上 


考え られ , Cluster 分 析 の 結果 が 複雑 に な っ て , 解 
釈 が 困難 に な っ た 場合 に 有効 で ある . 

各 被 験 者 6 人 に つい て 分 析 し た 結果 , 設定 され た 
2 つの グル ー プ 間 の イン ター バル は , 17msec か ら 
24msec に わた り 平 均 は 20.6msec CH), 平均 値 
に つい て は , 従来 の 研究 結果 や ! の や と ほぼ 一 致す る も 
の で あっ た . 従来 の 研究 に お いて は , 2 つの グル ー 
プ 間 の イン ター ベル の 条件 を , 平均 値 な ど に より, 
統一 され た 値 と し て 設定 し て いた が , 本 研究 に お い 
て は , 個々 人 に つい て 設定 値 を 求め た . この 方 法 に 
ko, G.D. の 分 析 に お いて , 個人 差 を 考慮 し た 客 
観 的 な 分 析 が 行え を る も の と 考え られ る . 

Fig. 6 は , Cattell23 が 示し た , 情報 分 析 に お 
ける 測定 値 の 基準 と その 方 向 を 表わし た 図 で ある . 
従来 の 研究 に お いて は , 個体 間 の 変動 を 用 い 個体 内 
の 要因 を 分 析 す る も の (R-technique, 横断 的 研 
究 ) , 1 個体 に つい て の 時 間 ご と の 変動 を 用 い 個 体 
内 の 要因 を 分 析 す る も の (P-technique, 縦断 的 研究 
と 呼ば れ 本 分 析 で 用 いた ) が 主 で あっ た . し か し , 
この よう な 2 次 元 的 な 統計 処理 か ら は , 各 要 因 の 状 
況 し か 見 る こと が 出来 ず , 解析 方 法 と し て , 不 十分 
な 面 が 残っ て いる と 言え る . 

今後 は , 単に 状況 を 見 る だ け の 解析 体 係 を 一 歩 進 
DC, 「 必 要 な 条件 を 再現 する 方 向 に 情報 と その 解 
釈 が 蓄積 され る こと が 必要 で ある ] と 言う 情報 処理 
の 考え 方 を 導入 する 必要 が 生じ る で あろ う . 今回 用 
いた Cluster 分 析 に こつ いて は , こう 言っ た 情報 処理 
の 考え 方 を 取り 入れ , 個体 間 x 時 係 列 x 個体 内 の 3 
次 元 の 情報 を 一 時 に 処理 する 3 Dimension Cluster 
分 析 の 開発 "や が 進め られ て お り , 今後 の 成果 が 期待 
され る 所 で ある . 

更に , 臨床 医学 や 基礎 医学 分 野 に お ける 筋 電 図 学 
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Fig. 6. Three dimension analysis of data matrix 
(Cattell. 1946). 
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的 研究 に お いて G.D. の 生理 学 的 意義 , と くに 体温 
調節 上 の 意味 づけ や その 形成 機 序 に 関し て は 幾多 の 
報告 O め め が あ る . し か る に 今回 の 実験 結果 の 如く , 
等 尺 性 筋 収縮 の 経過 中 に 誘発 され る 筋 疲労 の 筋 電 図 
波形 に も G.D. が 記録 きれ , か つ ,G.D. の 形成 に あ 
ず か る 連続 する 二 個 の 筋 放電 ズバ ルス 間隔 が 20 msec 
以下 で ある 事 を 客観 的 に 証明 し た 研究 報告 は 今後 の 
体温 調節 研究 分 野 に 新しい 指針 を 写 を える も の で あり , 
臨床 医学 へ の 応用 も 可能 で ある ・ 


<2 約 


常温 環境 下 に て 6 名 の 健常 被験 者 の LPI 
(M. Biceps brachii) に 随意 ・ 等 尺 性 筋 収縮 を 誘起 
し , その 経過 中 に 出現 する 筋 疲労 の EMG 波 形 を デ 
ー タ 処理 シス テム (ATAC450 及 び PC-8001) を 用 
いて 解析 し 以下 の 結果 を 得 た . 

① 等 尺 性 筋 取 縮 に ょ る 筋 疲労 時 の EMG 波形 に は 
寒冷 ふる え 時 の EMG と 同様 に 群 化 放電 (grouping 
discharge: G.D.) が 観察 され た . 

②G.D. の 形成 機 序 を 解明 する た め に , G.D. の 
構成 成分 で ある 連続 する 2 個 の 筋 放電 バル ス の イン 
ター バル に つい て 相関 行列 を 求め 更に P-mode 
Cluster 分 析 を 行っ た . その 結果 得 ら れ た dendro- 
gram iX, 2 個 の バル ス の イン ター バル が 1 一 19 
msec £20~30 msec の 2 群 に 分 か れ , DO, mW 
の ヒス ト グ ラ ム は 筋 疲労 の 進行 過程 に お いて 相反 的 
な 動態 を 呈 し た 

③② で 示 し た Cluster 分 析 の dendrogram で , 
イン ター ベル が 2 群 に 区 分 で き な か っ た 例 に 関し て , 
相関 行列 を 求め た 後に Direct Varimax 法 に よる 因 
子 分 析 を 行い 因子 行列 を 作成 し た . 

因子 分 析 の 結果 か ら , G.D. 形 成 の イン ター バル 
は 19 msec~21 msec が 適当 値 と 判明 し , この 事 は 従 
来 の 研究 結果 が 妥当 で ある 事 を 立証 し た 事 に な る ・ 

DHID, 従来 の 研究 で は G.D. の 条件 設定 に 主観 
的 要素 が 混入 する 可能 性 を 残し て いた が , 本 研究 (上 
述 ②③) の 結果 は , 多 変量 解析 の 駅 使 に よっ て G・D. 
の 条件 設定 が ょ り 客 観 的 と な っ た 点 が 進歩 で ある . 

⑤ 体 熱 産生 の 寒冷 ふる を と 随意 ・ 等 尽 性 筋 収縮 の 
誘発 法 や , その EMG 波 形 分 析 , 特に , Cluster 分 
析 や 因子 分 析 の 統計 学 上 の 意義 に つい て 二 ・ 三 の 検 
討 が 加え られ た . 
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